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ABSTRACT 

Interactions between a supermassive black hole binary and the surrounding accretion disc can 
both assist the binary inspiral and align the black hole spins to the disc angular momentum. 
While binary migration is due to angular-momentum transfer within the circumbinary disc, 
the spin-alignment process is driven by the mass accreting on to each black hole. Mass trans¬ 
fer between different disc components thus couples the inspiral and the alignment process 
together. Mass is expected to leak through the cavity cleared by the binary, and preferentially 
accretes on to the lighter (secondary) black hole which orbits closer to the disc edge. Low 
accretion rate on to the heavier (primary) black hole slows the alignment process down. We 
revisit the problem and develop a semi-analytical model to describe the coupling between 
gas-driven inspiral and spin alignment, finding that binaries with mass ratio q < 0.2 approach 
the gravitational-wave driven inspiral in differential misalignment: light secondaries prevent 
primaries from aligning. Binary black holes with misaligned primaries are ideal candidates for 
precession effects in the strong-gravity regime and may suffer from moderately large (^ 1500 
km/s) recoil velocities. 

Key words: accretion, accretion discs - black hole physics - galaxies: evolution - galaxies: 
nuclei 


1 INTRODUCTION 

Following a galaxy merger, the supermassive black holes (BHs) 
hosted ( [Kormendy & Richstone|1995^ by the two merging galax¬ 
ies sink towards the centre of the newly formed stellar environ¬ 
ment through dynamical friction, forming a binary jBegelman e t al.| 
|1980[[Mayer|2013| ). Binary BHs can merge if the astrophysical en¬ 
vironment provides a way to dissipate their angular momentum in 
less than a Hubble time. Scattering from stars can bring the binary 
only down to parsec scales jPrank & Rees|19'7^ , below which the 
available phase space is quickly depleted, thus stalling the inspi¬ 
ral process {final-parsec problem, |Milosavljevic & Merr itt||2001 [ 
|Yu|2002^ . While triaxiality in the stellar potential may help driving 
the inspiral in elliptical gas-poor galaxies ( [Merritt & Poon|2004[ 
[Berczik et al.|2006l l, the interaction with gaseous disc(s) may actu¬ 
ally solve the final-parsec problem in gas-rich galaxies jArmitage| 
|& Natarajan|2002[|Cuadra et al.|2009| . Indeed, dynamical friction 

against a gaseous background can promote the merger bringing the 
binary separation to distances of the order of 0.1 pc within a time- 
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scale of 10-50 Myrs ( jEscala et al.|2005[|Dotti et al.|2007||2009| (. 
Further shrinking of the binary can proceed through what is known 
as type II migration in the context of planet-disc interaction. How¬ 
ever, such disc-assisted migration can only be effective at separa¬ 
tions smaller than ~ 0.01 pc, beyond which the disc becomes self- 
gravitating and will likely fragment and form stars (jLodato et al.j 


2009 1 . Finally, if the binary reaches separations close to ~ 10 ^pc, 


gravitational-waves (GWs) quickly become an extremely efficient 
way to drive the binary to a prompt merger (jPeters & Mathewsj 
[19631 |Peters|T964| i. Asymmetric emission of GWs in the late in¬ 
spiral and merger may imprint recoil velocities to the remnant BHs 
( Redmount&ReesM9^ l which can be as high as 5000 km/s jC 


panelli et al.||2007[ [Gonzalez et al^|2007[ [Lousto & Zlochower 


2011 1 , possibly larger than the escape speed of the host galaxies 
( [Merritt et al.[2004[ . 


A variety of electromagnetic signatures has been proposed 
to detect supermassive BH binaries, which however remain elu¬ 
sive jPotti et al.|2012[[Schnittman|2013[[Bogdanovic|2015[ . The 
most convincing evidence comes from double active galactic nu¬ 
clei (AGN) imaging, with the notable example of the radio galaxy 
0402-(-379 showing two compact cores with estimated separation 
of 7.3 pc jRodriguez et al.|2006[ (. Indirect evidence for supermas¬ 
sive BH binaries at sub-parsec scales suffer from higher uncertain- 
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2 D. Gerosa et al. 


ties; we mention in particular the case of the blazar OJ287 jVal-| 
|tonen et al.||2008) , where ~ 12 yr periodic outbursts have been 
interpreted as signature of a BH binary orbital motion. Identify¬ 
ing recoiling BHs trough observations is even more challenging 


(|Komossa|2012 

1 , but tentative candidates are nonetheless present 

([Komossa et al. 

2008l[Civano et al.|2012|[Decarli et al.|2014|[Koss[ 


mergers are the main target of all future space-based GW observa¬ 
tories. The eLISA (evolved Laser Interferometer Space Antenna) 
mission, recently approved by the European Space Agency, is ex¬ 
pected to detect hundreds merging binaries per year casting new 
light on our understating of such systems jSeoane et al.|201^ |Ba-| 
[rausse et al.|2015^ . 

Hydrodynamical interactions can not only assist the binary in¬ 
spiral, but are responsible for the reorientation of the two BH spins. 
BH spin alignment has a crucial impact on the merger dynamics 
and on the cosmic growth history of supermassive BHs. Strong re¬ 
coil velocities can only be achieved if the merging BHs are highly 
spinning and the two spin vectors are strongly misaligned between 
each other. Highly recoiling BHs can be significantly displaced 
from the galactic nucleus or even ejected from it. This has strong 
consequences on the coevolution of BHs and their host galaxy, 
with mild recoil regulating the BH growth and large kicks veloci¬ 
ties strongly affecting the feedback process ( [Blecha & Loeb|2008| >. 
If large recoils causing BH ejections are present, this affects the 
fraction of galaxies hosting supermassive BHs jSchnittman|200 7l 
|Volonteri et al.|2070||Gerosa & Sesana|2015^ and consequently the 
predicted (e)LISA event rates i Sesanaj 2007| ). From the GW data 
analysis point of view, spin misalignments introduce a richer struc¬ 
ture in the expected signals that carries precious information on the 
binary dynamics (|Vecchio|2004||Klein et al.|2009||0’Shaughnessy| 
|et al.|20r^ and can improve the parameter estimation process by 
up to an order of magnitude ( [Chatziioannou et al.||2014) >. At the 
same time, since accurate waveform modelling is required in GW 
searches, spin precession makes the waveform generation more 
challenging ( [Pekowsky et al.|[2013[ [Klein et ^|2014| >, dramati¬ 
cally increasing the parameter-space dimensions that need to be 
explored. 

It is thus important for both electromagnetic and GW obser¬ 
vations to understand in which region of the BH-binary parameter 
space we expect significant spin misalignment, which is the goal 
of this paper. The physical process responsible for the reorienta¬ 
tion of BH spins during the long gas-driven inspiral is the so-called 
Bardeen-Petterson effect ( [Bardeen & Petterson||1975| l, where the 
general relativistic Lense-Thirring torque between the BH and a 
misaligned disc warps the accretion disc and secularly aligns the 
BH spin with the disc angular momentum. The Bardeen-Petterson 
effect does not just affect the binary dynamics during the gas-driven 
phase, but leaves a deep imprint in the subsequent GW-driven inspi¬ 
ral, where precession effects are strongly dependent on the residual 
misalignments left by the astrophysical environment ([Kesden et al.[ 
[2010a|b[ [Gerosa et al.|2013l l. It is therefore an essential ingredient 
to predict the spin configuration at merger. 

The effectiveness of the Bardeen-Petterson effect in aligning 
the spins to the binary angular momentum within the time-scale of 
the merger has been recently investigated by multiple authors. [Bog-[ 
[danovic et al.[ ( [2007t first made some order-of-magnitude estimates 
of the alignment time for a single BH with its own disc and found 
that it is much smaller than the merger time concluding that align¬ 
ment is likely in a gaseous environment. A similar conclusion was 
obtained by [Doth et al.[ ( (2010[ (, who found short alignment time- 
scales of ~ 2 X 10'’ yr. One notable achievement of the study per¬ 


formed by [Dotti et al.[ ( [20T0t is combination of smoothed-particle 
hydrodynamics (SPH) simulations and a semi-analytical treatment 
of the Bardeen-Petterson effect jPerego et al.|2009y through which 
they have been able to quantify the residual misalignment to either 
10° or 30° depending on the gas temperature. In a previous work 
( [Lodato & Gerosa|2013t , some of us revisited these estimates con¬ 
sidering the previously neglected effects of non-linear warps in the 
misalignment propagation trough the disc. Their conclusion is that 
the alignment time can be significantly longer than lO’ yr if the 
initial misalignment are large, thus casting doubts on the ability of 
the disc to align the binary. [Miller & Kjolik[ ( |20T^ made a sensible 
step forward, pointing out that the spin alignment process in BH 
binaries may actually be sensibly faster than for isolated BHs, be¬ 
cause of the stabilizing effect of the companion that increases the 
degree of disc warping close to the holes. 

In this paper, we argue that the binary mass ratio plays a key 
role in estimating the spin-alignment likelihood in merging BH bi¬ 
naries. We present a semi-analytical model to compute the inspiral 
and alignment processes from the properties of the circumbinary 
disc. On the one hand, the mass ratio strongly affects the binary in¬ 
spiral rate, marking the onset of different disc-morphology regimes 
when either a gap or a proper cavity can be opened. On the other 
hand, and perhaps most importantly, the binary mass ratio sets the 
amount of differential accretion on to the two components of the 
binary system. The Bardeen-Peterson effectiveness in aligning the 
spins depends sensitively on the mass accretion rate through each 
single disc. We quantify this quantity constructing prescriptions 
based on results of hydrodynamical simulations. Gas is expected to 
preferentially accrete on to the lighter binary member that therefore 
aligns faster. Accretion rates on to the heavier BH is consequently 
smaller and may prevent it from aligning. 

This paper is organized as follows. Our model to compute the 
relevant time-scales is introduced in Sec. [^ where both the inspi¬ 
ral and the alignment problems are treated. Sec. describes our 
main findings, namely the role of the binary mass ratio in the spin- 
alignment process, and its relevance for cosmologically motivated 
binary distributions and kick predictions. Finally, in Sec.j^we draw 
our conclusions and stress the possible caveats of our analysis. 


2 BINARY AND DISCS MODELLING 


We model the gas environment surrounding merging BH binaries 
through three different accretion discs: mass may reach the binary 
from galactic scales forming a circumbinary disc, and later be ac¬ 
creted on to the individual BHs from circumprimary and circum- 
secondary discs respectivel}^ We define R to be the binary sep¬ 
aration, Ml and M 2 to be the two BH masses (with Mi > M 2 ), 
Mbin = Ml +M 2 to be the total mass of the binary, q = M 2 /M 1 ^ 1 to 
be the binary mass ratio and S, = UjGMf jc to be the spin of any of 
tw o BHs (where i = 1,2 an d dimensionless spin 0 ^ a, ^ 1). Follow- 


Syer & Clarke (1995 1 , we also define a measure of the locan 
circumbinary-disc mass at a generic radius r as M(r) = 47rE(r)j^ 


mg 


' The name minidiscs can also be found in the literature to indicate circum¬ 
primary and circumsecondary discs. 

^ For typical values of the surface density exponent, this is a reasonable es¬ 
timate of the rigorous value obtained radially integrating the surface density 
jShakura & Sunyaev|1973t . In addition, the migration rates are set by the 
properties of the disc in the vicinity of the planet, and the local nature of 
this parameter is thus more relevant than the total disc mass. 
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Spin alignment in BH binaries 3 


and finally qi\sc(r) = M(r)/Mbm to be the disc-to-binary mass ra¬ 
tio. When studying the spin alignment, we will refer to the mass 
of the aligning BH as M (meaning either Mi or M 2 ) and the mass 
of the other binary member as Me. While the accretion rate of the 
circumbinary disc Mbin determines the inspiral process, the align¬ 
ment time-scales are only determined by the rates Mi and M 2 at 
which mass reaches the circumprimary and circumsecondary discs 
respectively. It is also useful to define / to be the dimensionless 
value of Mbin in terms of the Eddington accretion rate 

j.Mbln 

Mbin — f — , (1) 

fEdd 

where fedd = Knc/AirG ~ 4.5 x 10*yr is the 
and Ke is the opacity for Thomson electron scattering. 

In this section, we first present a new estimate for the inspiral 
time-scale by interpolating estimates computed in different regimes 
and we discuss the circumbinary disc self-gravity condition to eval¬ 
uate such inspiral time-scale. Secondly, we summarize the main 
findings of |Lodato & Gerosa| ( |2013) on the spin-alignment time- 
scale and we explore the effect of the companion on the individual- 
discs structure. We finally model mass transfer and differential ac¬ 
cretion on to the different discs. 


Salpeter (1964 1 time. 


2.1 Gas-driven inspiral 

If BH mergers do happen in nature, it is likely that the gas-driven 
phase is the bottleneck of the whole binary evolutionary track. 
Therefore, the time spent by the binary in such phase gives us an 
estimate of the total time available to align the BH spins through en¬ 
vironmental processes before merger. Although gas-driven inspiral 
is mediated by the torques exerted by the disc on to the binary, a de¬ 
tailed description of the torques is not necessary to correctly derive 
the migration rates. Ultimately, the migration rate is controlled by 
the rate at which the disc is able to redistribute the angular momen¬ 
tum gained from the binary, and the torques will adjust to give the 
correct rate (e.g.,|Armitag e|20T0) . This same mechanism is called 
type II migration in the context of protoplanetary discs (|Lin & Pa-| 
|paloizou|1986[ ). Depending on the ratio between the BHs and the 
circumbinary disc masses, we identify three possible regimes: 


HSyer & Clarke||1995[ [Ivanov et al.||T^ |Lodato et alT][200^ 
[Baruteau & Masset|2013^ and reads 


M2^Mm) 

M{R) 


(3) 


which correctly reduces to tv{R) in the limit M 2 ^ M{R). 

(iii) For comparable mass binaries (Mi < M 2 ), the secondary- 
BH potential cannot be neglected. The gap at the secondary loca¬ 
tion now becomes a proper cavity in the disc with radius ~ 2R 
cleared by both BHs fMacFadyen & Milosavljevic|2008| ). A rich 
phenomenology may be present, including disc asymmetry and 
growing eccentricity, and can only be captured using hydrodynam- 
ical simulations ( [Cuadra et al. ||2009[ |Roedig et al. pOfTj |2012t 
| Shi et al.|2012^ . An approximate expression for the inspiral time- 
scale in the comparable mass regime has been presented by |Rafikov| 
l |2013^ , assuming the binary potential to be represented by a New¬ 
tonian potential produced by the binary total mass. He obtains 


tin = 


M1M2 
Mbin M(R) 


tu(R) , 


(4) 


where the correction factor Mi /Mbin models the expected speed up 
due to the higher angular momentum flux induced by the binary 
mass. The same mass-ratio dependence has been very recently ob¬ 
tained by |Dotti et al.| ( (20751 1 integrating the torque at the edge of a 
2R-wide cavity 

Here, we propose a smooth analytical interpolation between 
the time-scales obtained in the different regimes given by 


tin = 


Ml M2+M(R) 
^ M(R) 


UR), 


(5) 


which correctly reduces to either Eq. l[^, or in the relevant 
limits. Various numerical factors in Eq. Q -as already acknowl¬ 
edged by IRafiko^pO 13 [ l himself- and different possible definitions 
of the viscous time-scale may modify our this estimate of a factor 
~few. The accretion rate of the circumbinary disc Mbin only enters 
in the merger time-scale through the viscous time-scale, which can 
be rewritten as 


UR) = 


3 M(R) 

4 ’ 


( 6 ) 


(i) For small mass ratios (M 2 ^ M(R) <C Mi), the secondary BH 
perturbs the disc of the primary, which reacts opening a gajj^at the 
binary separation ( |Lin & Papaloizou|1979t[Artymowicz & Lubow| 
|1994| >. Tidal interactions between gas particles and the secondary 
BH transfer angular momentum to the disc, thus decreasing the bi¬ 
nary separation. The secondary BH behaves like a fluid element in 
the disc, evolving at the viscous rate ( |Armitage & Natarajan|2002) 

^2 

tin=t„(R)~ - , (2) 

V 

where v = aCsH is kinematic viscosity coefficient of the disc, usu¬ 
ally jShakura & Sunyaev|1973^ rescaled to a dimensionless coeffi¬ 
cient a with the speed of sound c., and the disc height H. 

(ii) If the secondary BH mass becomes comparable to the disc 
mass (M 2 ~ M(R) Mi), the disc cannot efficiently redistribute 
the momentum acquired from the binary. The shrinking rate con¬ 
sequently decreases. An analytical expression for the inspiral time- 
scale in this regime can be computed directly from the angular- 
momentum conservation equation in the thin-disc approximation 

^ Binaries with very small mass ratios (g < 10^, see |Armitage|20l6| can¬ 
not open a gap; however, such low mass ratios are not expected to be rele¬ 
vant in the supermassive BH context. 


since both M(R) = 47rR^E(R) and Mbin — 37rizS(R) are related to the 
surface density E of the circumbinary disc in the stationary limit. 
Combining Eqs. Q, § and 1^, we obtain our final estimate of 
the inspiral time-scale, to be compared to the individual alignment 
time-scales. 


tin = 


3 (tq) qdisciR) q tedd 


(l+?)2 


/ 


(7) 


While the low mass ratio regime is relatively well tested in the 
planetary community (e.g. [Nelson et al.|2000[|Bate et al.|20d3] but 
see below for possible caveats), the regime of high mass ratios has 
not been explored extensively. At the moment only few simulations 
of disc driven migration of a binary have been conducted (|Mj^| 

[Fadyen & Milosavljevic|2008[|Cuadra et al.|2009[|Shi et al.|2012| >, 
which test only a small part of the parameter space. For example, 
for a ratio q=l /3, [Cuadra et al.| ( (2b09^ find in 3D SPH simulations 
a migration rate of R/R = -2x 10“^ Q, where Q = (GMun/R^Y^^ is 
the orbital frequency of the binary. They compare this value with 
the analytical formula from [Ivanov et al.[jl999^ , ou r Eg. ([^, which 


yields R/R = -3 x 10 ^ fl, in very good agreement. MacFadyen & 


[Milosavljevic[ ( [2008) find in 2D grid-based simulations of g = 1 bi¬ 
naries that the inspiral time-scale is roughly the viscous time-scale, 
reduced by gdisc, which is consistent with Eq. 0. The simulations 
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described so far neglect the details of the angular momentum redis¬ 
tribution mechanism, which in the standard picture is the magneto- 
rotational instability (MRI,| ^lbus & Hawley|1991^ , and typically 
adopt the Shakura & Sunyae^ 1973| > a parametrization (e.g., |Ma^ 
|Fadyen & Milosav ljevic|2 008t in o rder to reduce the computational 
cost. Only recently |Shi et al.| ( |2012[ l were able to perform global nu¬ 
merical simulations of migrating binaries including the MRI. They 
found that magnetohydrodynamics effects slightly enhance the mi¬ 
gration rate with respect to the purely hydrodynamical case (a fac¬ 
tor of ~ 3 when compared to |MacFadyen & Milosavljevic|2008| >. 
They also observed that the accretion of material with a higher spe¬ 
cific angular momentum than the binary can make the binary gain 
angular momentum, which however is offset by the higher torques 
they measure from the disc. Given the number of other uncertain¬ 
ties present in the model, we are thus satisfied that our expressions 
can be used reliably. 

Recent numerical simulations jPuffell et al.|2014{|Durmann| 
|& Kley|2015^ in the planetary community have questioned the va¬ 
lidity of type II migration, casting doubts that a regime where the 
satellite behaves like a test particle exists at all. In particular, the 
simulations show that it is possible to achieve faster (up to a fac¬ 
tor of 5) migration rates than what expected from Type II theory. 
These simulations have only been run for a fraction of a viscous 
time, and it is still unclear if this result holds on the time-scale of 
the merger. For this reason we neglect these results in what follows, 
and note that this makes our estimates an upper limit for the merger 
time-scales. 

Finally, we note that the simulations conducted so far, to the 
best of our knowledge, have explored relatively thick discs, with an 
aspect ratio ranging from 0.05 to 0.1 (e.g.,|MacFadyen & Milosavl-| 
[jevic|2008|[Cuadra et al.|2009HRoedig et al.|2012[|Shi et al.|2012f 

This is significantly thicker than the value we derive in the next 
section and it is not clear how the results would change with more 
realistic values (cf. Sec.|^. 


2.2 Self-gravity condition 


The inspiral time-scale reported in Eq. ([^ depends on the binary 
separation R. For typical disc structures ^Shakura & Sunyae v|1973 [ 
|Goodman||2003|>, tm is a steep monotonically increasing function 
of R ( |Haiman et al.||2009) . Most time will be spent by the bi¬ 
nary at large separations, while the remaining inspiral is completed 
rather quickly. The time available to align the spins -which the 
spin-alignment time must be compared to- is roughly the inspiral 
time-scale tin evaluated at the largest separation of the disc-driven 
evolution. 

A natural physical limit on the size of the circumbinary disc 
is set by the disc self gravity. Local gravitational stability under 
axisymmetric disturbances is guaranteed up to the fragmentation 
radius Rf where the|Toomr^s (|1964|l parameter equals unity: 


Q = 


Csil 

ttGE 


( 8 ) 


At separation R > Rf, self gravity cannot be neglected and the disc 
is gravitationally unstable (cf. e.g. |Lodato|2007| for a review). The 
evolution of gravitationally unstable discs has been investigated in 
great details in recent years ( [Lodato & Rice|2004[|Rice et al. 12005 ) 
ICossins et al.|[200^ . If the cooling time is smaller or of the or¬ 
der of the dynamical time ( |Gammie|2001| >, the disc will fragment 
into gravitationally bound clumps, although the actual fragmenta¬ 
tion threshold is debated ( |Meru & Bate|20f^ . For values appro¬ 
priate to AGN discs, the disc is expected to fragment, create stars 


and thus deplete the area surrounding the binary of gaseous ma¬ 
terial, possibly halting the inspiral jLodato et al.|200^ . Using the 
vertical-equilibrium equation Cj/U = //, the self-gravity stability 
condition Q = 1 can be rewritten as 


^disc(R) ■ 


M(R) 

4/bin 



(9) 


evaluated at R = Rf. The fragmentation radius is likely to lie in the 
outer region of the circumbinary disc, dominated by gas-pressure 
and electron-scattering opacity ( [Shakura & Sunyaev|1973^ . Assum¬ 
ing viscosity to be proportional to the gas pressure (/3-disc) and 
setting the mass of the accreting object to Mbin, one gets for the 
fragmentation radius [ Shakura & Sunyaev| 1973[ |Goodman |2003[ 
[Haiman et al.|2009) 


Rf ~ 10 


~0.05 


5 GMbin f Mbir 


c2 V 10''4fo 
f Mb.„ 

V lO’Mg 







pc. 


( 10 ) 


For a separation r in such region, the disc aspect ratio reads 



In this paper, we evaluate the inspiral time-scale of Eq. 0 at 
the fragmentation radius: R = Rf. This is a rather conservative as¬ 
sumption, being tm monotonically increasing with R ( [Haiman et al.j 
[2009[ l and being Rf the largest separation at which gas can be found 
under the form of a circumbinary disc. We are therefore assuming 
-somehow overcoming the final parsec problem- that some previ¬ 
ous mechanisms are efficient enough to shrink the binary separation 
down to Rf. 

From Eqs. izini we find that the inspiral time tin scales only 
mildly with the viscosity a and the binary total mass Mbin. As for 
the accretion rate /, the implicit dependence from Eq. fu} and ( [Tol l 
is also mild; the explicit dependence 1// in Eq. 0 is still present 
but will cancel when compared to the spin-aligment time (Sec. in}. 
On the other hand, the dependence on q plays a crucial role when 
comparing each spin-alignment time with the inspiral time-scale, 
and is therefore the main subject of this study. 


2.3 BH spin alignment 


The circumprimary and cirbumsecondary discs interact with the 
BHs through the Bardeen-Petterson effect. [Bardeen & Petterson[ 
( [1975[ showed that a viscous disc initially misaligned with the 
equatorial plane of a spinning BH naturally relaxes to a coplanar 
state in the inner regions, while the outer disc may retain its origi¬ 
nal misalignment. [Ree's |( |7?78) realized that, by Newton’s third law, 
the outer disc must react by pulling the BH towards complete align¬ 
ment (or antialignment) of the spin with the orbital angular momen¬ 
tum of the outer disc itself. 

Angular momentum is initially transferred from the spin to 
the inner disc trough relativistic Lense-Thirring precession and fi¬ 
nally to the outer disc by the propagation of warps, i.e. vertical 
shearing by close, misaligned, gas rings (Scheuer & Peiler|1996[ 


[Lodato & Pringle|2006[ [Martin et al.|20ci7y r Warp propagation i 

ruled by a vertical viscosity coefficient V 2 , which is generally differ¬ 
ent than the kinematic viscosity coefficient v introduced above. As 
done for a, let us introduce a vertical-viscosity coefficient Q 2 such 
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that V 2 = aiCsH ( [Papaloizou & Pringle ||1 9 83| >. In the small-warp 
-which in our case actually means small-misalignment- limit, the 
warp-propagation coefficient is related to the kinematic viscosity 
by jPringle|1992[|Ogilvie|199^ 


0.2 


( 12 ) 


1 4(1-H7a^) 

2a 4 + a^ ’ 

and, in particular, it is independent of the misalignment ip between 
the inner disc and the outer disc. A full non-linear theory of warp 
propagation has been computed by Ogilvie| ( |l999| > and later veri¬ 
fied numerically by |Lodato & Price|j[2010^ . Non-linearities intro¬ 
duce a qualitatively new dependency on ip, which can lower the 
value of a 2 by a factor of ~ 7 for large misalignment angles (see 
Fig. 1 in|Lodato & Gerosa|2013f. In this paper we consider the full 
non-linear expression a 2 (a, (p) as derived by |Ogilvt^ ( |1999^ , which 
reduces to Eq. © for < 1 . 

Lense-Thirring precession efficiently aligns the disc up to the 
Bardeen-Petterson radius Rbf, defined to be the disc location where 
the inverse of the Lense-Thirring precession frequency ( [Wilkins | 

[jwnj 


Q.LT(r)-2 


equals the warp propagation time = 
2/3 /H \ -4/3 


^Iv 2 , i.e 




02 


(13) 


(14) 


For a single BH-disc system, Rbp coincides with the maximum 
warp location (warp radius) Rw and marks the boundary between 
the (quickly aligned) inner disc and the (still aligning) outer disc. 
The time-scale over which the outer disc finally aligns the BH 
spin can be found by computing the torque acting on the disc at 
Rw ^Natarajan & Pringle|1998| l. A single BH of mass M and di¬ 
mensionless spin a aligns with the angular momentum of the sur¬ 
rounding accretion disc within (jScheuer & Feiler|1996]|Natarajan| 


|& Pringle|1998[|Lodato & Gerosa|2013 l 



M f 

a H 

4i ~ 3.4 

— a 

— 


M V 

02 r 

where 

M is 

the 


2/3 


(15) 


mary/circumsecondary disc and H/r its aspect ratio evaluated at 
the warp radius. We note here that the alignment time 41 is sensibly 
smaller that the growth time M/M for reasonable viscosities 
a ~ 0.1 and aspect ratios H/r ^ 0.001. BH mass and spin 
magnitude can be therefore considered fixed during the alignment 
process ( [King & KolbTT^ see Sec.|^. 

The Bardeen-Petterson effect can drive the BH spin towards 
either alignment or antialignment with the outer disc. [King et ah] 
([2005^ showed that the system antialigns if 


6»>7r/2 and L(Rw)<2S, 


(16) 


where 9 is the angle between the BH spin and the angular momen¬ 
tum of the outer disc, L(Rw) is the angular momentum of the inner 
disc (i.e. inside the warp radius) and S is the BH spin. The BH spin 
aligns with the outer disc if any of the two conditions above is not 
satisfied. Once 9 is provided (cf. Sec. [^, the misalignment p be¬ 
tween the inner-disc angular momentum and the outer-disc angular 


4 The warp-propagation coefficient actually depends on the radial deriva¬ 
tive of the local inclination of the disc ip, see Eq. (1) injLodato & Pricej 
Here we implement the same approximation ip ^ pas already done 
by Lodato & Gerosa|^2013^. 


momentum is given hy p = 9 m the aligned case, while p = n-9 
if the system tends towards antialignment. Note that even in cases 
where the BH spin antialigns with its own disc, the net effect is al¬ 
ways to reduce the misalignment with the binary plane (|King et al.| 
[20^ . 

2.4 Effect of the companion on disc-spin alignment 

So far we have only considered the alignment of a single BH 
with its surrounding accretion disc. Here we discuss the effect of 
a far (R Rw) companion on the alignment process. Such effect - 
neglected in our previous study 1 Lodato & Gerosa|2013[ (- has been 
recently pointed out by [Miller^^^^oir^ 201 ^ in the supermas- 
sive BH binary case, while Martin et al.j ( 2009^ have previously 
considered the same interaction for stellar-mass BHs with stellar 
companions. For a further study, see jTremaine & Davis| ( |2014l (. 

If the aligning BH is part of a binary system, the gravitational 
potential felt by an orbiting gas ring is perturbed by the presence of 
the companion ( jLubow & Ogilvie|2000[ [Ogilvie & Dubus|2001| (. 
The binary gravitational potential can be expanded in series of r/R 
(where r is the distance of the gas ring from the BH and R is the 
binary separation; see e.g. jKatz et al.|198^ : to leading order, the 
resulting torque is perpendicular to the angular momentum of the 
gas ring L, causing precession about the angular momentum of the 
binary Lbin- The precessional frequency can be obtained by averag- 
ing the torque over the binary orbital period and read^ jPetterson 

wn\ 

4 I?3 \gm) 


(17) 


where /3 = |Lbin • L|, M is the mass of the aligning BH and Me is the 
mass of the companion. Note that in our notation Me = qM in the 
case of the primary BH, but Me = M/q when the alignment of the 
secondary is considered. 

If a spinning BH is part of a binary system, both Lense- 
Thirring and companion-induced precession are present. The com¬ 
panion drives the system towards alignment with the angular mo¬ 
mentum of the binary, which tracks the plane of the circumbinary 
disc (see Sec. on this point). At the same time, the inner disc 
is being aligned to the BH spin by Lense-Thirring precession. In 
practice, the companion reduces the frame-dragging efficiency: ma¬ 
terial could stay misaligned with the BH spin at closer locations, 
thus speeding the alignment process up jMiller & Krolik||2013| (. 
This effect can be quantified by computing the locations at which 
the two contributions are equally important. The Lense-Thirring 
time Qp/ equals the warp propagation time at Rbp, as given by 
Eq. (g. On the other hand, the disc is now expected to be max¬ 
imally warped at the warp radius Rw, where the Lense-Thirring 


contribution matches the companion one Qlt = He 1 Martin et al. 


2009 ) 1 . From Eq. (fTS) and (|T7) one gets ^Martin et al. 
& Krolik|201^ 


2009 Miller 


8a M\ 


^'^=\WWe 


i?2/3 


(f)"’ 


(18) 


If Rw sj Rbp, the companion term can be neglected and the closer 
location at which misaligned material can be found is still ~ Rbp- 
The alignment speed-up discussed by [Miller & Krolikj p013[ l is 

3 We do not quote the sign of the precession frequency, because it only sets 
the precession direction about Lyn which is not important for our order-of- 
magnitude estimate. 
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6 D. Gerosa et al. 



Figure 1. Effect of the companion on the location of the wai'p radius. Primary and secondary BHs are considered on the left-hand and right-hand panel 
respectively. Contours map the ratio between the maximum-warp location (where the external torque from the companion matches the Lense-Thirring 
torque) and the Bardeen-Petterson radius Rgp (where the Lense-Thirring time equals the warp propagation time) derived for an isolated BH. The dependence 
on the binary mass ratio q and the spin of the aligning BH (either or ^ 2 ) is reported on the x and y axes. A companion BH is expected to speed up the 
alignment process up to a factor of ~ 2 in most of the parameter space, meaning a 0(1) uncertainty in the estimate of the alignment time. This figures have 
been produced taking Mbin = IO^Mq, H/r = 0.001, / = 0.1, a = 1 /2a.2 = 0.2, /3 = 1 and R = Rf [see Eq. l |l9[ ]. 


relevant if Rw ^ Rbp, because warped regions are present closer to 
the hole. AtR = Ri (cf. Sec. |2.2[ t, we find 

io’Mq J 

HA y'" I-JAV”'" 

(0.001/ VO- 1 / vO. 2 / \l/ 2 aj ■ 


i<BP 


,- 2 / 9 ^- 4/9 


M+Mc 


2/3 


Fig .[^ shows the dependences of Rw/Rbp on the binary mass ratio 
and the spin magnitude of the aligning BH, both for primaries and 
secondaries. Slowly rotating BHs are less affected by the presence 
of a companion because the spin set the magnitude of the frame¬ 
dragging term. For fixed total mass Mbin = M + Mc, primaries are 
more sensible to the companion than secondaries, because their 
gravitational radius is larger and Lense-Thirring precession can be 
matched more easily by the additional precession term. 

In this paper we use the simple expression reported in Eq. ( |15^ 
to compute the spin-alignment time, as formally obtained for an 
isolated BH-disc system. Our analysis [Eq. jl9| l and Fig. 0 shows 
that the position of the warp radius can be modified by a factor 
of ~ 2 if th e BH is part of a binary system. The alignment time 


fai oc iNatarajan & Pringle|l998 


Miller & Krolik 


2013 I can 


therefore only be lowered by a factor of ~few. From Eq. 
this assumption may not be valid if (i) the binary is very massive 
Mbin -Sj IO^Mq, (ii) the individual discs are thinner than the cir- 
cumbinary disc at the fragmentation radius H/r < 0 . 001 , (iii) the 
binary accretion rate is close to the Eddington limit / > 0.1. A 
more complete understating of the alignment process in BH binary 
systems requires explicit integrations of the angular momentum 
equation ^Martin et al.|200^ . This goes beyond the scope of this 
work, which instead focuses on getting an estimate for the align¬ 
ment time-scale. 


2.5 Cavity pile-up and differential accretion 

The accretions rates of the individual discs Mi and M 2 depend on 
the circumbinary-disc accretion rate Mbin, since the formers are fed 
by the latter. Here we develop a simple prescription to link these 
three quantities. 

Accretion from the outermost regions of the circumbinary disc 
on to the binary BH is suppressed because of either the perturba¬ 
tion of the secondary (low-i? regime) or the two-body central po¬ 
tential (high-^ regime). Therefore, the binary may only accrete at a 
lower accretion rate Mi +M 2 ^ Mbm. Mass tends to pile up at the 
outer edge of the cavity created by the binary itself: accretion -and 
therefore spin alignment- is still possible if gas streams can pene¬ 
trate the cavity and reach the BHs. We call Mgap the mass accretion 
rate that overcomes the cavity pile-up: this gas will sooner or later 
accrete on to either the primary or the secondary BH, i.e. 

Mgap = Ml + M 2 . (20) 


Gas-stream propagation is an intrinsically multi-dimensional 
non-linear phenomenon that requires dedicated hydrodynamical 
simulations to be studied in detail. In particular, the dynamics of 
gas accretion through the cavity is strongly dependent on the binary 
mass ratio q, since qualitatively different regimes are present. |Mac-| 
|Fadyen & Milosavljevic| ([2008|l first discussed equal-mass binary 


et al. 


simulations, while the ^ = 1/2 case has been presented by Hayasaki 
1 2007|l and extended to 17 = 1/3 by Cuadra et al. 1 2009^ and 


Roedig et al. 1 2012 1 . A recent major improvement has been made 


by D’Orazio et al. 1 2013^ and |Earris et al.| ( [2014| > who extensively 
studied the dependence on q of the mass rate overcoming the cavity 
pile-up. 

[D’Orazio et al.| ( |2013) present 2D hydrodynamical simula¬ 
tions in the range 0.003 5/ ^ ^ 1 assuming fiducial values a = 0.01 
and ///r = 0.1. Accretion on to the binary is indeed limited to nar¬ 
row gas streams and it is typically suppressed by a factor of 2-5 
when compared to a single-BH disc of the same mass. They detect 
the presence of two physical regimes: 
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Spin alignment in BH binaries 1 




Figure 2. Numerical fits to hydrodynamical simulations to compute the accretion rates of the two members of a BH binary. Left-hand panel shows the fraction 
of the accretion rate Mgap/Mbin fhat penetrates trough the disc cavity and reaches one of the two binary members. A quadratic interpolation is performed to 
numerical results by |D’Orazio et al.]^2013) here reported in Eq. (D- On the right-hand panel, the binary accretion rates is broken down between the primary 
and the secondary BHs. The ratio of the accretion times in the systems simulated by|Farris et al.H2014J appears to be described by the simple prescription 
(M 2 /M 2 )/(Mi/Mi) = Ijcp in Eq. { 22 } . 


(i) For high mass ratios 0.05 < ^ ^ 1, the presence of the binary 
strongly modulates the streams. Streams are generated by devia¬ 
tions from spherical symmetry in the binary potential: more asym¬ 
metry is present for equal mass binaries that therefore show less 
mass pile-up at the cavity edge and more binary accretion. In such 
regime, the ratio Mgap/Mw„ is expected to increase with q. 

(ii) In the low-mass ratio regime q < 0.05, the secondary BH 
quickly swipes through the disc accreting most of the material com¬ 
ing from large distances: a single gas stream is present feeding the 
secondary BH. Such effect gets more pronounced when the mass 
ratio is lower and consequently the ratio Mgap/A/bin decreases with 
q. 


Such results are shown in the left-hand panel of Fig. |2 were the 
“MidAr-LoA(^” simulations by |D’Orazio et al.| ( |2013 1 are con¬ 
sidered. The minimum in Mgap/Mbin separates the two physical 
regimes just described above. We interpolate the results from the 
simulations performed by |D’Orazio et al.|j2013^ with the ansatz 

=po +pi log(^) -H p 2 log^ (q) (21) 

Mbm 


and best-fitting coefficients po = 0.8054, pi = 0.9840 and p 2 = 
0.3818. For a given circumbinary-disc accretion rate Mbin, Eq. lH^ 
specifies the mass rate Mgap which overcomes the cavity pile-up and 
accretes on to either the primary or the secondary BH. The subse¬ 
quent study by [Farris et al.|pOT4l l found that the ratio between the 
total accretion rate on to either one of the two BHs and the accre¬ 
tion rate on to a single BH of the same total mass may exceed unity, 
thus casting doubts on whether such fraction can be interpreted as 
Mgap/Mbin- To bracket this uncertainty, we use Eq. 1 21 1 as our ref¬ 
erence model but we al so s tudy an additional variation where we 
fix Mgap = Mbin (cf. Sec. 


,__ M- 

[Farris et al.j ( [2014^ recently performed 2D grid simulations 
(assuming H/R = Q.\ and q = 0.1), specifically addressing the feed¬ 


ing of the individual discs from streams penetrating the cavity. 
They systematically find that the secondary BH accretes faster than 
the primary, mainly because the former orbits closer to the cav¬ 
ity e^e. Their results are here reported in the right-hand panel of 
Fig. M where the ratio of the accretion times M,/M, (t = 1,2) is 
showed as a function of the binary mass ratio q. Symmetry implies 
Ml ~ M 2 for binaries with high mass ratios, while lower values of q 
show pronounced differential accretion in favour of the secondary. 
A qualitatively different regime is detected for the lowest of their 
simulated cases q = 0.025: the cavity is not efficiently cleared by 
the secondary BH, and mass from the circumbinary disc directly 
flows inwards forming a large circumprimary disc. As pointed out 
in Sec. [2.1[ such change in the dynamics of the system is expected 
for lower mass ratio, where the disc should form a small annular 
gap rather than a large hollow cavity. To directly reach the circum¬ 
primary disc, gas should be able to flow past the secondary escap¬ 
ing its gravitational attraction. As pointed out by [Farris et al.[j20T4l l 
themselves, the actual turning point in q is likely to be highly de¬ 
pendent on the thickness of the disc and possibly on the viscosity. 
As shown recently by [Young et al.[ ( [2015^ in the context of binary 
stars, direct flowing from the circumbinary to the circumprimary 
disc is easier for thicker discs, where the stronger pressure forces 
can make part of the material “skirt” the Roche lobe of the sec¬ 
ondary, eventually reaching the primary Roche lobe and being cap¬ 
tured by its gravitational attraction. Due to such uncertainties, in 
this work we deliberately ignore the onset of such low-g regime 
when considering differential accretion. The growth-time ratio pre¬ 
sented b ^Farris et al.[ j2014[ l appear to be well approximated by 
(see Fig.|^ right-hand panel) 


M2/M2 _ 1 
Ml/Ml ~ q- ' 


( 22 ) 


Due to such pronounced differential accretion, the prescription 
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Figure 3. Combined effect of mass pile-up at the edge of the disc cavity 
and differential accretion in unequal-mass binaries. Prescriptions for the 
accretion rates presented in Eqs. j2H23) and Fig. are here summarized. 
As a function of q, we show the fraction of the accretion rate of the outer 
circumbinary disc Mbin that (i) overcomes the cavity-pile up and accretes 
on to the binary (Mgap, solid); (ii) is captured by the secondary BH which is 
clearing the cavity (M 2 , dashed); (iii) finally accretes on to the primary BH 
(Ml, dotted). 


here presented may formally predict super-Eddington rates for the 
secondary BHs in the low-g regime. This has no relevant im¬ 
pact on our model: defined the secondary Eddington ratio to be 
fi = etEddAf 27^2 = (l+q)M 2 tf /qMhm, the Eddington limit /2 = 1 
is only marginally reached for very high circumbinary-disc accre¬ 
tion rates (/ ~ 1) and low mass ratio q < 0.05 (assuming a typical 
accretion efficiency e ~ 0.1). Lower values of / shift the critical 
mass ratio at which the Eddington limit is formally reached to even 
lower values. Finally, we note that, as in Sec. o the thickness 
values explored by the simulations considered in this section are 
significantly higher than those expected for massive-BH binaries 
(Sec.[g. 

From Eq. ( |20^ and ^22\ one gets 

Ml = , M 2 = j:|:^"gap . (23) 

Fig. combines the results presented in Eq. l |21[ ) and Eq. j23[ l. 
For equal-mass binaries, ~ 80% of the incoming mass may accrete 
on to the binary and it is equally distributed between the two bi¬ 
nary members. Unequal-mass binaries present differential accre¬ 
tion Ml < M 2 that grows stronger as q is decreased. When q is in¬ 
creased from q ~ 0.2 to unity, gas streams start to flows towards the 
primary and M 2 /Mhm consequently flattens. On the other hand, if q 
is decreased from q ~ 0.05 to 0.001, the secondary orbit gets closer 
to the inner edge of the cavity jD ’Orazio et al.|201^ : more mass 
can overcome the cavity pile-up and it is almost entirely accreted 
by the secondary. 


3 RESULTS: DIFFERENTIAL MISALIGNMENT 

In this section we compare the spin-alignment time and the inspi¬ 
ral time. We first outline the regions of the parameter space where 


Parameter 

Fiducial model 

Synthetic distributions 

<1 

Free parameter 

Power-law distributions 

^bin 

Not relevant 

Not relevant 

fll, 02 

1 

Either 1 (E) or 0.1 (C) 

9l, 02 

Extremize over 

Random variables 

R 

Rf (fragmentation) 

Rf (fragmentation) 

H/r 

0.001 

0.001 

OL 

0.2 

0.2 

f 

Not relevant 

Not relevant 

Ml , M 2 

Eqs. j2l| and j23| 

Eqs. l|2l| and |23| 


Table 1. Choice of the binary and disc parameters in our time-scale compar¬ 
ison for both the fiducial case (Sec. ED and the cosmologically motivated 
distributions (Sec. |3.2t . 


misalignments are foreseen (Sec. ED; secondly, we fold our model 
into synthetic supermassive-BH binary populations (Sec. |3.2| l; and 
we finally present a preliminary study to address the impact of our 
findings on the occurrence of large post-merger kicks (Sec.|3.3[l. 


3.1 Misaligned primary BHs 

3.1.1 Fiducial values of the parameters 

The circumbinary disc properties enter the inspiral time An while 
the primary/secondary alignment times tai are set by individual-disc 
parameters. The ratio 4i/fin in general depends on the binary sep¬ 
aration R, the three disc aspect ratios H/r, the gas viscosity a, the 
accretion rates of the circumbinary Mwn and the individual discs 
M1.2, the BH masses Mi and M2 (or equivalently q and Mbin), the 
orientation angles 9i and 62 , and the BH spin magnitudes ai and fl 2 - 
We first specify a fiducial model by taking likely values of all these 
parameters and we later perform a small parameter study around 
such model. Table 13 . iTTI summarizes the values we assume for the 
parameters, highlighting the differences with the next section. We 
discuss our choices as follows. 


• As detailed in Sec. |2.2[ a rather conservative assumption can 
be made by evaluating the inspiral time at the fragmentation radius 
Rf. This is a measurement of largest separation where the inspiral 
can be driven by interaction with a gaseous environment and is typ¬ 
ically believed to be the bottleneck of the whole binary evolution. 

• In our fiducial model we fix the aspect ratios of all discs to 
H/r = 0.001. As reported in Eq. 0 for the circumbinary disc, 
the aspect-ratio dependences on the other parameters (namely the 
viscosity, the accreting mass and the accretion rate) are not crucial 
to evaluate the inspiral time-scale, and will be here neglected for 
simplicity. For the same reason, we assume the individual discs to 
share the same aspect ratio of the circumbinary disc (cf. the analo¬ 
gous assumption made by |Miller & Krolik|2013^ . 

• Unless specified, we fix a = 0.2. A parametric study on the 
viscosity has already been presented previously ( [Lodato & Gerosa] 
|2013| > and the alignment process has been found to be overall quite 
independent of a. Negative azimuthal viscosities are formally pre¬ 
dicted by the non-linear warp propagation theory for a < 0.1 and 
large misalignments ip ( |Ogilvie| 1999 [[Lodato & Gerosa|20f3) : the 
evolution of the disc in these cases is unclear and out of the scope 
of this study (see [Nixon & King][2012[ and [Nixon et al.|2013] for 
extensive discussions). 

• As described in Sec. 2^ the BH accretion rates Mi and M 2 are 
related to the circumbinary disc accretion rate Mbin, conveniently 
expressed through the dimensionless quantity / in Eq. Q- In our 
fiducial model we implement Eqs. ED and ( |23| l. Once H jr is fixed, 
the alignment likelihood Li/fin is independent of / because both 
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Figure 4. Comparison between the alignment times and the inspiral time 
as a function of the binary mass ratio q. Gas interactions have time to align 
the BH spins with the orbital angular momentum of the binai^ if the align¬ 
ment time /ai (shaded areas) is smaller than inspiral time t\a (dashed line). 
While secondary BHs (lower, green area) aligns for all mass ratios, this 
is not the case for the primary members (upper, red area) which retain 
their initial misalignments if q is low enough. Our fiducial model is as¬ 
sumed here: H/r = 0.001, a = 0.2, / = 0.1, and maximally spinning BHs 
= 1- Warp non-linear propagation theory introduces uncertainties 
(thus the shaded areas) of a factor of ~ 2.5 in the alignment times. The 
con'esponding times obtained with the linear theory are shown with dotted 
lines for comparison and always underestimate the non-linear result. 


times scale as 1 // [cf. Eqs. and a . This is a point of improve¬ 
ment over our previous estimate ( [Lodato & Gerosa|2013| l, where an 
effective dependence on / was introduced when decoupling the in¬ 
spiral and the alignment processes. For concreteness, the overall 
scale of Fig.|^below is computed assuming / = 0.1. 

• Within our assumptions, both the inspiral and the alignment 
times are independent of the binary total mass Mbm. This is com¬ 
patible with |Haiman et ai:] ( |2009| l when is evaluated at the frag¬ 
mentation radius. 

• The orientations angles 6 i and 62 set the warp efficiency 02 
(Sec.|2.3| and their dependence is the main point raised by |Lodato| 
|& Geros^ ( |2013[ l. In the following, we bracket such uncertainties 
extremizing 4i over all possible orientations. 

• For simplicity, we consider maximally spinning BHs 
(fli = 02 = 1) unless specified. The status of supermassive-BH spin 
measurements has been recently reviewed by [Reynold^ ( |2013| l: 
some highly spinning BHs are found, but the current statistic is too 
low to provide a complete picture of the spin magnitude distribu¬ 
tions. The effect of the spin magnitude on the alignment likelihood 
can however be easily predicted, because the alignment time scales 
as fai oc cf. Eq. 1 15 1. 


3.1.2 Predicted time-scales 

The key dependence of the problem is the one on the binary mass 
ratio q, which both marks the onset of different inspiral regimes and 
sets the importance of differential accretion. Fig.|^shows the inspi¬ 
ral and the alignments times as a function of q for our fiducial set of 
parameters. The uncertainty in the initial misalignments 9i causes 
the alignment times to appear as stripes in the figure, rather than 


Variation 

Primary BH 

Secondary BH 

Fiducial 

[0.14-0.23] 

(0.14) 

[n-n] 

(n) 

0 = 0.3 

[0.17-0.28] 

(0.18) 

[n-n] 

(n) 

Q = 0.4 

[0.19-0.32] 

(0.21) 

[n-n] 

(n) 

0 = 0.5 

[0.21 -0.35] 

(0.23) 

[n-n] 

(n) 

H/r= 

[0.07-0.12] 

(0.07) 

[n-n] 

(n) 

Hjr^ 10^2 

[0.28 - 0.48] 

(0.29) 

[n-n] 

(n) 

H/r= 10“' 

[0.60- 1] 

(0.61) 

[0.02 - n] 

(n) 

a, = 0.2 

[0.09-0.14] 

(0.09) 

[n-n] 

(n) 

a, = 0.5 

[0.12-0.19] 

(0.12) 

[n-n] 

(n) 

11 

p 

CO 

[0.14-0.22] 

(0.14) 

[n-n] 

(n) 

^bin ~ ^gap 

[0.07-0.12] 

(0.07) 

[n-n] 

(n) 


Table 2. Binary mass ratios marking the transition between aligned and 
misaligned spins. For any variat ion from our fiducial model [a,- = I, a = 0.2, 
H/r=\ 0“^, Mgap given by Eq. j2l| ] , we report values q such that BH spins 
in binaries with q < q are expected to be left misaligned (i.e. 4i > fj^) by 
gaseous interactions. Values in square brackets refer to the lower and upper 
limit of q due to the initial-misalignment uncertainty foreseen using non¬ 
linear warp propagation. Values in round brackets show the analogous result 
when the linear theory is considered, and notably underestimates the value 
of q. Misaligned secondaries are typically not present (as indicated with 
“n”) unless some of the parameters are cranked up to unrealistic values. 


lines. For comparison, we also show (dotted lines) the behaviour 
predicted by the linear warp-propagation of Eq. © where 41 is 
independent of Oj. The linear theory underestimates the alignment 
time by up to a factor of ~ 2.5, as already pointed out by |Lodato| 
|& Geros^ j2013| l. Fig. illustrates the main result of this paper: 
while secondaries are found aligned (t-^.i <C fin) for every value of 
q, primary BHs only align if ^ > 0.2. Light secondaries may pre¬ 
vent primaries from aligning. If such BHs were misaligned before 
the disc interactions, these misalignments are carried over to the 
next stages of the binary evolution. As briefly explored in Sec. |3.3| 
this differential alignment between the two binary members will af¬ 
fect the subsequent GW-driven inspiral, the merger phase and the 
properties of the remnant BHs opening for the possibilities of large 
kicks. 

A short parametric study around our fiducial model is shown 
in Table|^ where we compute the values of q which mark the onset 
of the misaligned regime (i.e. where 4i = fin). As expected ( [Lodato | 
|& Geros^|2013| (, the alignment process is rather independent on 
a with thresholds varying from ^ ~ 0.17 to 0.35 if a is increased 
from 0.2 to 0.5. Notably, the alignment likelihood is also rather in¬ 
dependent on the spin magnitudes ai and < 22 , because of the mild 
scaling of 4i [cf. Eq. (|15H. Alignments times are longer for maxi¬ 
mally spinning BHs ai = a 2 = 1 chosen for our fiducial model, but 
misaligned primaries are predicted for mass ratios fir ~ 0.15 even 
when moderately spinning BHs are considered. Perhaps more sur¬ 
prisingly, the alignment process appear to be strongly dependent 
on the discs aspect ratio H/r which enters linearly in fin and with a 
lower power in 4i. Only primaries with q > 0.6 have enough time to 
align their spins in thicker discs ///r ~ 0.1, even when maximally 
rotating BHs are considered. Moreover, if H/r is large enough, the 
inspiral time may become comparable with the secondary align¬ 
ment time within the physical uncertainty due to initial spin orien¬ 
tation. As already pointed out (Sec. se also Sec.|^below), the 
disc thickness is one of the main uncertainties in the current mod¬ 
elling of binary-disc interactions. Details of the gas streams leaking 
through the disc cavity have also a notable effect: the largest value 
of q where misalignment is foreseen drops down to ~ 0.12 if all of 
the gas of the circumbinary disc ends up being accreted by either 
one of the two BHs (i.e. if Mbin = 


■ Mgap, cf. 


Farris et al. 


20141. 
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Figure 5. Mass ratio distributions in the synthetic supermassive-BH binary 
populations developed by |Arun et al.H200^ . Four models are available, for 
different prescriptions of the accretion geometry (Efficient versus Chaotic) 
and the BH seeds (Large versus Small). Data sets are binned within the 
range q S [0.01,1] and fitted with power laws N oc q'^. The best-fitting 
spectral indexes 7 are reported in the legend for each model. The histogram 
normalization has been inherited from the original models and is irrelevant 
to our purposes. 

3.2 Cosmogically-motivated distributions 

Our findings are relevant if supermassive-BH binaries with spins 
and mass ratios in the misaligned regime are present in nature 
and detectable. While electromagnetic observations already con¬ 
strained almost a hundred supermassive BH masses jMcConneli] 
|& Ma|[MT3) and a handful BH spins ( |Reynolds||2013| l, the mea- 
surements of the global properties of the supermassive-BH binary 
population is the main goal of future space-based GW observa¬ 
tories. eLISA ( [Seoane et al.pofst will detect hundreds binaries 
per year up to redshift z ^ 10 with sufficient signal-to-noise ratio 
0(10- 100) to measure accurately both individual-source parame¬ 
ters and their statistical distribution jBarausse et al.|20T5) . 

Here we present a simplified analysis to address whether the 
misaligned-spin regime highlighted above is relevant in such con¬ 
text. Publicly availably synthetic distributions of merging BH bi¬ 
naries have been developed by the LISA collaboration in the con¬ 
text of the LISA Parameter Estimation Taskforce ( |Arun et al.|2009^ 
and later updated by |Sesana et al.|f2011^ . The authors developed 
four merger-tree models of BH evolution, varying over only two 
ingredients, considered to be the main sources of uncertainty. 

(i) The mass of the BH seeds. In the small seed scenario (S), first 
BHs of mass ~ lOOM© are initialized as remnants of Population 
III stars at z ~ 20 and evolved according to |Volonteri et al.|p003| l. 
In the large seed scenario (L), BH with mass IO^Mq are formed 
from gaseous protogalactic discs at z ~ 15 to ~ 10 as developed by 
|Begelman et al.| ( [2006| l (see also |Lodato & Natarajan|2006^ 

(ii) The accretion geometry. If accretion efficiently (E) occurs 

® http://www2.iap.fr/users/volonter/LISA_catalogs.html 



Figure 6. Fraction of BH spins in binary systems that align within a factor r 
of the inspiral time as predicted using the publicly available distributions 
by |Atun et al.|j2009^ . Alignment predictions using the model presented in 
this paper can be read at r = 1 : all four distributions show that ~ 8 % of 
BH primaries may fail to align during the gas-driven inspiral, while strong 
differential accretion quickly aligns all secondaries. Fractions P at larger 
and lower values of r predict the alignment likelihood in case of systematic 
modelling errors on either the inspiral or the alignment time. 


on few long episodes, the BHs will generally be spun up during 
their cosmic evolution ( |Thorne|1974| l. On the other hand, accretion 
may also happen to be chaotic (C), on many short episodes (|King| 
|& Pringle|2006^ . In this case, lumps of material accreted in random 
directions spin on average the holes down. 

This approach results into four models, referred to as SE, SC, LE 
and LC. Eig.j^shows the extracted mass ratio distributions, together 
with power-law fits A oc in the range q € [0.01,1]. We obtain 
7 = 0.17 for LC, 7 = 0.56 for LE, 7 = 0.40 for SC and 7 = 0.63 
for SE. The spin-magnitude distributions presented by |Arun et ^ 
( |2009^ are strongly peaked towards slowly spinning BHs for the C 
models and a ~ 1 for the E models. This is a direct consequence 
of their simplified accretion treatment, which is either completely 
coherent or completely chaotic; broader distributions are predicted 
for more realistic evolutionary models where such assumption is 
relaxed ( |Barausse|2012[[Dotti et al.|2013[[^sana et al.|2014^ . Spin 
orientations are not tracked during the cosmic evolution by |Arun| 
|et al.|p()09| l: spins are assumed to efficiently align in models E, 
while their directions are kept isotropic in models C. 

Fig.|^show the cumulative fraction of aligned BH using these 
four synthetic BH-binary population^ We sample the mass ratio 
q over the fitted power-law distributions from Fig. spin magni¬ 
tudes are set to a\ = 02 =0.1 in the C models and a\ = 02 = \ in 
the E models, to mimick the strongly peaked distributions of|Arun| 
|et al.lpOO^ . For simplicity, we fix the disc properties to our fidu¬ 
cial values [H/r = 0.001, a = 0.2, Mgap given by Eq. l |2l| ] and we 
sample over a uniform distribution in cos 6 i to extract values of the 
alignment time within the initial-orientation uncertainty presented 

^ We are aware of the inconsistency of our procedure, being the binary 
mass ratio distributions used here coupled to the spin orientations: at each 
merger tree level, the properties of the daughter BHs do depend on the spin 
orientations of their progenitors [cf. Sec. 2.1 in [Gerosa & $esan^j2015) 
and references therein]. 
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in Sec. ED Fig.j^shows, for each value of r, the fraction of binaries 
P for which 4i < AnT. Sections at r = 1 correspond to our current 
model: while all secondaries align during the inspiral, up to ~ 8% 
of the primaries may not have time to align their spin before merger. 
This statement appear to be rather independent of the population 
synthesis model chosen. In particular we find P(4i,i < An) = 0.93 
for LC, P(h\,\ < fin) = 0.96 for SC, P(fai.i < An) = 0.92 for LE and 
P(fai.i < All) = 0.93 for SE. Two main effects are here combined: 
while the E models present higher spin magnitudes (hence longer 
alignment times) than the C models, they also predict a steeper pro¬ 
file in the mass ratio (Fig. |^, with fewer small-^ binaries (hence, 
on average, shorter alignment times). Those misaligned, rapidly ro¬ 
tating, primaries BHs predicted by the E models are ideal targets 
for strong-gravity precession effects in the late inspiral and merger 
( |Schnittman|2004[|Kesden et al.|2015] Sec. |3.3[ >. 

Fig.j^also provides intuitions on the consequences of system¬ 
atic errors in our time-scale estimates. If the inspiral (alignment) 
time is larger (smaller) of a factor r = 10, all binaries in the sample 
align by the end of the gas-driven inspiral. On the other hand, if the 
inspiral (alignment) process is 10 times faster (slower), i.e. r = 0.1, 
only 60%-80% of the primaries aligns. 


3.3 Differential misalignment and kick velocity 


The most notable consequence of our findings is a clear prediction 
for the spin-orientation angles at the onset of the GW-driven in¬ 
spirals: a non-negligible fraction of supermassive BH binaries ap¬ 
proaches the GW-driven phase with 9i^Q and 62 — 0. 

If the binary lies in the same plane of the circumbinary disc 
Pvanov et al.|1999[ [Miller & Krolik|2013[ see Sec. 0’ the angles 
9i may be taken as estimates of the misalignment between the BH 
spins and the binary angular momentum Lbin, and used to estimate 
the properties of the post-merger BH. While final mass ( [Barauss^ 
|et al.|20i^ and spin ( [Barausse & Rezzolla|2009} do not critically 
depend on the spin misalignments, these are crucial to predict the 
final recoil ( [Campanelli et al.]|2007[ [Gonzalez et al.pOOV^ . The 
largest kick velocities (up to ~ 5000 Km s“‘) are attained for maxi¬ 
mally spinning, equal-mass BH mergers with moderately large mis¬ 
alignments 9i ~ 50° ( jLousto & Zlochower|201 1||2013^ . 

Here we perform a preliminary study to estimate the impact 
of our findings on the kick velocity distribution. To maximize the 
effect, we consider maximally spinning BHs in binaries with mass 
ratio q = 0.2, right at the onset of the misaligned regime highlighted 
in Sec. [33 (cf. Fig.|4|. Numerical-relativity fitting formulae are 
available to compute kick velocities, but precession effects during 
the GW-driven inspiral must be taken into account, especially for 
configurations with sensibly different spin tilts 9i ^ 92 jSchnittman| 


2004||Kesden et al.|2010a|b||Kesden et al.|2015||Berti et al.|2012[ 


Gerosa et al. |2013| l. GWs start driving the merger at the decou- 
plin g radius ((Armitage & Natarajan||2002[ [Gold et aL]|2014[ see 

alsojMilosavljevic & Phinney|2005) 


rdtc — 760—^— — 


(S 


-2/5 


H/r 

0.001 ) (1-h#/5 ’ 


(24) 


where the angular momentum losses in GWs dominate over the 
viscous evolution of the disc. We first transfer the spin orientations 
from the initial separation rdec = 760GMbi,i/c^ to rgn = 10GMbin/c“ 
using the precession-averaged formalism recently presented by 



cos 9i at rdoc 


Figure 7. Maximum kicks velocity vyck (colour scale and marker size) as a 
function of the misalignment angles cos 0 , = Si ■ Lbin, measured at the decou¬ 
pling radius rjec (x and y axes for primary and secondary BHs respectively). 
We consider maximally spinning BH binaries with q = 0.2, predicted to be 
at the onset of the misaligned regime unveiled by our astrophysical model. 
Large kicks are foreseen in the 61 y 0, 02 — 0 region, predicted to be as- 
trophysically relevant. Superkicks with vj. > 2000 km/s are not likely in 
gas-rich environments because binaries with larger mass ratio are expected 
to align before mergers. 


Kesden et al.| ( |2015 1 ; and we finally appljQtte numerical-relativity 


fitting formula by Lousto & Zlochov(^ l |2013l l at r^n- 

Fig. prelates the spin orientation at the decoupling radius to 
the maximum kicks velocities allowed in each configuration. No¬ 
tably, higher kicks are found in the 71 ' 0 , 02 — 0 region, which 
we predict to be populated by the Bardeen-Petterson effect (~ 8 % 
of the cases from the models used in Sec. |3.2| l. The evolution of 
the spin orientation in the GW-driven inspiral can be qualitatively 
understood in terms of two special families of configurations {spin- 
orbit resonances, \Sc\mMmdin\lQQ^ , in which the projections of the 
two spins on the orbital plane are either aligned or anti-aligned to 
each other. Kicks are suppressed (enhanced) by spin-precession ef¬ 
fects for binaries in the aligned (anti-aligned) family ( [Kesden et al.j 
[2010b[ (. Configurations lying in the 0i 7 / 0, 02 — 0 region of the 
parameter space are likely to be attracted into the anti-aligned fam¬ 
ily (cf. e.g. Fig. 5 in |Gerosa et al.|2013^ : as binaries approach the 
merger phase, most of their spin-precession cycles is spent with the 
two spins forming an angle ~ tt when projected on to the orbital 
plane. These configurations are qualitatively similar to the standard 
superkick configuration ( [Campanelli et al.]|2007[ [Gonzalez et al.| 
[2007[ l and notably predict high kick velocities. 

Our prediction is that, if merging rapidly rotating BHs are 
present in gas-rich environments, kicks as large as vt ~ 1500 km/s 
can happen. Such kicks can make the BH wander in the galaxy 


* We refer the reader to Sec. 2.1 in [Gerosa & Sesan^[2015l for details 
on the fitting formulae implementation. In their notation, we assume a ran¬ 
dom initial phase Ad) at rjec. In order to disentangle the dependence of the 
kick velocity on the spin orientations, we maximize over the orbital phase 
at merger 0 (thus only showing the maximum kick allowed in each config¬ 
uration). The relevance of 0 on the results presented in this section can be 
easily predicted because the kick velocity scales roughly as vj. oc cos 0 [c.f. 
e.g. Eq. (2) in jCampanelli et ^[2007[ ]. The kick velocity is independent 
of A/bin - 
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outskirt for times as long as 10-100 Myr with displacements of 
lO^^pc jOualandris & Merritt 2008 Komossa & Merritt 2008 Si 


[jacki et al.|201 l[|Gerosa & Sesana|2015[), possibly at the level of 

observational consequences ( Komossa|2012^ . Larger values of v* 
are only possible in merging binaries with mass ratio closer to the 
equal-mass case (for more quantitative information see e.g. Fig. 3 
in |Lousto & Zlochower|20f3) . Both BFls in these binaries are pre¬ 
dicted to be found aligned at merger ~ 02 ~ 0 (Sec. |3.l| l, which 
limits the kick velocity to ~ 300 km/s. Our analysis shows that 
superkicks with vr > 2000 km/s are disfavoured in gas-rich envi¬ 
ronments where the Bardeen-Petterson effect comes into play. 


4 DISCUSSION AND CONCLUSIONS 

Alignment of BH spins in merging BFl binaries may be differen¬ 
tial. Using a semi-analytical model, we find that light secondaries 
may accrete almost all mass leaking tough the binary gap and pre¬ 
vent primary BHs from alignment. In particular, such differential 
alignment occurs for binary with mass ratio q < 0.2. Gaseous inter¬ 
actions have enough time to align both spins in binaries with mass 
ratio closer to the equal-mass case (Sec. in}. We implement our 
analysis trough a time-scale argument, comparing the time needed 
to align the BH spins in the Bardeen-Petterson effect 4 i to the to¬ 
tal time available in the gas-driven inspiral phase The align¬ 
ment and the inspiral processes are coupled by the accretion rates: 
while the binary migration is set by the circumbinary disc mass 
rate, alignment is powered by the mass accreting on to each BH. 
Mass from the circumbinary disc is expected to pile up at the outer 
edge of the cleared cavity, suppressing the alignment process. On 
top of this, mass leaking trough the cavity is found to preferen¬ 
tially accrete on to the secondary BH which orbits closer to the 
disc edge. This causes the alignment time of the primary BH to be 
several orders of magnitudes longer than that of the secondary, and 
possibly even longer than the inspiral time. Differential accretion 
is a key, previously neglected, feature to tackle the spin-alignment 
problem: for comparison, |Miller & Krolik| ( |2013| > only quoted a fac¬ 
tor of ~ g between the alignment times of the two BHs. While 
powerful for its simplicity, our time-scale argument fails to cap¬ 
ture the dynamics of the alignment process: more elaborate models 
involving numerical simulations are needed to predict the resid¬ 
ual misalignment of primary BHs that cannot be aligned with the 
Bardeen-Petterson effect, and to estimate how close to complete 
alignment secondaries can be found in realistic environments. 

We present preliminary results to address the relevance of our 
findings on to the supermassive-BH cosmic history. Using publicly 
available synthetic populations, we find that binaries in differen¬ 
tial misalignment are expected in realistic cosmological scenarios 
(Sec. |3.2| l. A fraction of ~ 8 % of the BH primaries are found mis¬ 
aligned at merger even in models predicting large spin magnitudes, 
opening for the possibility of large kick velocities. Merging BHs 
with spin angles 02 ~ 0 and 0 i 7 ^ 0 are subject to the largest kicks 
velocities available for their mass ratio and spin magnitudes. In par¬ 
ticular, misaligned primaries in BH binaries with q ~ 0.2 may suf¬ 
fer kicks as large as ~ 1500 km/s, while higher mass ratios are 
needed to succeed in achieving proper superkicks (Sec. |3.3^ . Bi¬ 
naries approaching the merger phase with differentially misaligned 
spins will exhibit pronounced precession effects in the later GW- 
driven inspiral phase jSchnittman|2004l |Kesden et al.|2015| l. Or¬ 
bital plane precession modulates the amplitude of the GW cycles, 
encoding information of the astrophysical environment on to the 
emitted GW pattern ( [Gerosa et al.|2014| [Vitale et al.|2014[ Trifiro 


et al. in preparation). These features may in principle be used to 
recall and constrain our models using future space-based GW ob¬ 
servations, although more work is needed to quantify these state¬ 
ments. 

Several assumptions have been made in developing our mod¬ 
els, some of them worth of future improvements. First and perhaps 
most importantly, our model estimates whether spin misalignments, 
if present, are carried over towards merger. The same dynamical 
processes that bring the binary together may play a role in deter¬ 
mining the spin directions before Type II migration takes place. 
While star scattering is unlikely to affect the spin orientations be¬ 
cause it does not present any preferred direction, this may not be 
the case for previous larger-scale gas interactions. Dynamical fric¬ 
tion against gaseous environment may be crucial to promote the 
binary formation: even if short (~ 10 Myr;[Escala e t al.[ 2005[ [Dot ti] 
[ et al.[20071 l, this phase presents interesting dynamics involving tidal 
shocks and nuclear cusp disruption jVan Wassenhove et al.|20i4l l 
whose possible consequences on the spin directions still need to 
be explored. Secondly, our model only estimates whether the BH 
spins aligns to the angular momentum of the disc, while strong- 
gravity effects in latest inspiral and merger phase depend on the 
misalignments between the BH spins and the binary angular mo¬ 
mentum. The further assumption of alignment between the binary 
orbital plane and the circumbinary disc ( [Ivanov et al.[1999l ( is nec¬ 
essary to estimate the properties of the post-merger BHs, in par¬ 
ticular the kick velocity. Thirdly, we have neglected the BH mass 
growth during the alignment process. Differential accretion brings 
binary towards larger mass ratios on time-scales ~ MjM. While 
this effect can be safely neglected on the time-scale of the align¬ 
ment process 4 i ~ IQT^M/M [cf. Eq. 1 15 1 ], it may not be negligi¬ 
ble on the time-scale of the inspiral. However, this point may only 
be important for aligned binaries which do not present large kick 
velocities anyway. As extensively discussed is Sec. [2.4[ we also 
neglect the presence of the companion when estimating the align¬ 
ment time (cf. [Miller & Krolik|2013| l. However, Eig.[T] shows that 
this effect (a factor of ~ 2 in 4 i) mostly affects q^ \ binaries where 
both BHs aligns anyway. This point is worth further investigation, 
but sensible modelling efforts are likely to be required because the 
presence of two external torques (Lense-Thirring precession and 
the companion) cannot be fully captured within a time-scale argu¬ 
ment ( [Martin et al.|2009) . Einally, we have assumed that all values 
of q are allowed on cosmological grounds. Assuming the BH mass 
correlates with the galaxy mass, galaxy pairs with g < 0.1 may fail 
in forming close binaries because of strong tidal interactions before 
the galactic merger jTaffoni et al.|200^[Callegari et al.|2009[[Van] 


[Wassenhove et al.|2014 1 . Mass stripped away from the secondary 

galaxy may sensibly increase the delay time between the galaxy 
and the BH mergers, possibly even preventing the BH binary for¬ 
mation. 

We stress that the impact of the disc aspect ratio on the 
Bardeen-Petterson effect is still not understood and can potentially 
be crucial. Both migration process (Sec. 113 and gas streaming 
through the binary cavity (Sec. [2.5^ have only been simulated with 
sensibly thicker disc (typically H/r ~ 0.01 - 0 . 1 ) than those pre¬ 
dicted for discs surrounding supermassive BHs {H jr ~ 0.001, see 
Sec. [2.2[ l. In particular, a sensibly lower amount of gas may be able 
to leak trough the cavity in thinner discs, possibly slowing down 
the alignment process. Although we are aware of the computational 
constraints in simulating thin discs, we stress that such simulations 
are needed to validate the analytical expressions assumed here, and 
we point towards the importance of pushing these numerical efforts 
to lower values of the aspect ratio. 
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Merging supermassive BH binaries are unique systems where 
gravity and astrophysics both play together to shape the dynamics. 
BH spin alignment (or misalignment) is an imprint of angular mo¬ 
mentum transfer between the astrophysical and the relativistic side 
of BH binaries whose potential still need to be fully uncovered. 
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